Observation of universal behaviour of ultracold quantum critical gases 
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Quantum critical matter has already been studied in many systems, including cold atomic gases. 
We report the observation of a universal behaviour of ultracold quantum critical Bose gases in a 
one-dimensional optical lattice. In the quantum critical region above the Berezinskii-Kosterlitz- 
Thouless transition, the relative phase fluctuations between neighboring subcondensates and spatial 
phase fluctuations of quasi-2D subcondensates coexist. We study the density probability distribution 
function when both these two phase fluctuations are considered. A universal exponential density 
probability distribution is demonstrated experimentally, which agrees well with a simple theoretical 
model by considering these two phase fluctuations. 
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The nature of quantum criticality [lH3| driven by quan- 
tum fluctuations is still a great puzzle, despite of the 
remarkable advances in heavy-fermion metals and rare- 
earth-based intermetallic compounds, etc [j|. New un- 
derstanding of quantum criticality is widely believed to 
be a key to resolving open questions in metal-insulator 
transitions [5|, high temperature superconductivity ^] 
and novel material design, etc. Cold atoms in optical lat- 
tices provide a unique chance to not only simulate other 
strongly correlated systems |7| , but also study some mod- 
els unaccessible in solid state systems, particularly for the 
Bose-Hubbard model |8l4lO| . Despite of its complexity, a 
strongly correlated system in quantum critical regime is 
expected to exhibit a universal behaviour described by a 
certain physical quantity. 

Here we report the observation of universal behaviour 
for ultracold quantum critical Bose gases in a one- 
dimensional optical lattice. Density probability distri- 
butions of the released gases are measured for differ- 
ent depths of the lattice potential. It was found that 
the density probability follows a simple exponential law 
when the Bose gases reach the quantum critical region 
above the Berezinskii-Kosterlitz-Thouless (BKT) transi- 
tion 11-141. This universal behaviour can be well un- 



derstood in terms of our theoretical model considering 
both the relative phase fluctuations of quasi-2D subcon- 
densates and spatial phase fluctuations of individual sub- 
condensates above the BKT transition. The method of 
density probability distribution should provide a unique 
tool for identifying certain quantum phases of optical lat- 
tice systems. 

Ultracold Bose atoms in ID, 2D and 3D optical lattices 
are widely studied by the Bose- Hubbard model [^]. At 
zero temperature, there is a continuous quantum phase 
transition from superfluid (SF) to Mott insulator (MI). 
Because of the strongly correlated quantum behaviour, 
the quantum critical point (QCP) at absolute zero tem- 
perature distorts strongly the structure of the phase dia- 



gram at finite temperatures, leading to the emergence 
of an unconventional 'V-shaped' quantum critical re- 
gion 0, [l5l - [l9{ (Fig. 1). In analogy with a black hole, 
the crossover to quantum critical gas involves crossing 
a 'material event horizon', which implies strongly that 
the quantum critical gas has a simple and universal be- 
haviour [3|. 

In the commonly used 3D optical lattice systems, the 
atoms in a lattice site are strongly confined in all direc- 
tions, with a mean occupation number per lattice site 
of about 1 ~ 3 [lOJ. Thus the spatial phase fluctua- 
tions for the atoms in a single lattice site can be omit- 
ted, and the relevant theoretical calculations based on the 
Bose-Hubbard model [8|,|9| and Wannier function [20| can 
give quantitative descriptions of almost all experimental 
phenomena [7|. In contrast, for ultracold Bose atoms in 
a ID optical lattice, unique characteristics may arise due 
to the following two reasons: 

(1) Quasi-2D Bose gas: In an experiment of ID 
optical lattice system as ours, there can be hundreds of 
atoms in a lattice site. The Bose gas in a lattice site 
becomes quasi- two-dimensional (quasi-2D) |2lM23|. if the 
local trapping frequency uj^ of a lattice site in the lattice 
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direction satisfles the condition hwz » ksT 

(2) BKT transition: For such a quasi-2D Bose gas, 
the critical temperature in the occupied lattice site is 
given by T2D ~ ^± (^;C (2))^^^ /ks, with A^, being the 
atomic number in a lattice site. Well below T2D, the 
whole system becomes a chain of subcondensates. If there 
is no correlation between the subcondensates in different 
lattice sites, a quasi-2D gas undergoes a BKT transition 
at Tbkt = 72£)/4 [111, [13]. Beyond this critical value 
of Tbkt, it is favorable to create vortices in quasi-2D 
subcondensates, and the unbinding of bound vortices will 
lead to strong spatial phase fluctuations within each sub- 
condensate. 

Since the BKT transition is crucial to revealing the 
property of cold atoms in a ID optical lattice, its critical 



1.0 



0.8 



0.6 



0.4 



0.2 



\ critical temperature of quasi-2D Bose gas 
critical temperature of the combined potential 



Type-I QC region / 

^''" a^ universal exponential I 

a5 y [a_J (jgnsity probability distribution / 

\ / 

^ / Type-I 



\ 



BKT transition / 



Ml region 



SF region 



^ Type-I I QC ^ 

^ j«gion y Type-ll U\ region 



80 



60 



40 



20 



(a1) s=0 



-2 2 QCP4 6 8 

log,„(U/J) 

FIG. 1: Finite temperature phase diagram of ultracold Bose 
gases in a one-dimensional optical lattice. C/ and J repre- 
sent the strengths of the on-site repulsion and of the nearest 
neighbor hopping in the Bose-Hubbard model, respectively. 
By increasing the strength of an optical lattice, U j J can be 
increased by several orders of magnitude. The existence of 
the quantum critical point (QCP) at zero temperature leads 
to a 'V-shaped' quantum critical (QC) region. Apart from the 
ordinary SF-QC-MI transition, the BKT transition for quasi- 
2D subcondensates can also occur at sufficiently low temper- 
atures. In the QC region above the BKT transition (Type-I 
QC region), there exists a universal behaviour in the density 
probability distribution. All the phase boundary lines, except 
for the dashed line between QC and MI regions, are calculated 
with a total atom number of 1.1 x 10"" and trap parameters 
used in the experiment. Thar is the critical temperature for 
the purely harmonic trap in absence of the optical lattice. The 
gray line displays a phase transition path from SF to Type-I 
QC region. Data points of (a4)-(a6) along the gray line are 
obtained from the corresponding absorption images in Fig. 2 



temperature is also plotted in the phase diagram (red 
line in Fig. 1). The BKT transition line divides the 
QC region into two parts (Type I and Type II). Our 
experiments focus on the transition from SF region to 
QC region above the BKT transition. 

The experiments started with pure ^^Rb condensates 
confined in a magnetic trap with axial and transverse 
trapping frequencies of {w^,Wz} = 27r {83.7, 7.6} Hz. 
The ID optical lattice was formed by a retroreflected 
laser beam of A = 800 nm along the axis (z direction) 
of the condensate. This laser beam was ramped up to a 
given intensity over a time of 50 ms, yielding a lattice po- 
tential Vopt = sEr sin^ (27rz/A), with En being the recoil 
energy of an atom absorbing one lattice photon. After 
a holding time of 10 ms, we suddenly switched off the 
combined potential and allowed the cold atomic cloud to 
expand freely for a time of 30 ms. 

The expanded atomic gases were probed using the con- 
ventional absorption imaging technique. Since the probe 
beam is applied along x direction, what an absorption im- 
age records is the two-dimensional (y—z) column density 
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FIG. 2: Density probability distribution for different lattice 
depths. (al)-(a6) in the left column are absorption images 
showing the density distribution of the released atomic clouds. 
The field of view is 0.9 x 0.9 mm^, and the pixel size is 
A^ — 9.0 X 9.0 fira^. Tc denotes the critical temperature 
of the atomic gas in the combined trap. For moderate s, 
the interference fringes are prominent. Further increasing s, 
however, we see complete disappearance of the interference 
fringes (a6). This is due to the crossover from SF to QC re- 
gion as s is increased. (bl)-(b6) in the right column are the 
corresponding density probability distributions. Open circles 
are the data points calculated from the density distributions 
of the images in the left column. The blue solid curves give 
the exponential density probability distribution Pq defined 
in text. When s becomes large enough, the data points agree 
well with the exponential curve (see b5 and b6), which shows 
clearly a universal behaviour in the QC region above the BKT 
transition. 



profile of the atomic cloud. We denote by Nf {y, z) and 
-/Vf ' (y, z) the number of photons detected in the pixel 
at position (y, z) with and without the atomic cloud, 
respectively. Then the density distribution is written 
as n2D (y, z) = In [iVf (y, z) /Nf (y, z)] \^/a. 
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FIG. 3: Density fluctuations of the released atomic clouds. 
(a)-(d) are one-dimensional cuts through the corresponding 
images. For s = 30.7, the density distribution is highly fluctu- 
ated along both the horizontal (z) and vertical (y) directions. 



being the absorption cross section of a single atom and 
A the pixel size. 

The left column of Fig. 2 ((al)-(a6)) displays the den- 
sity distributions for different lattice depths s. The tem- 
peratures of the system were inferred from the conden- 
sate fractions as described in the Appendix. From the 
top three images ((al)-(a3)), we see an increasing of the 
side peaks in the interference patterns as s is increased. 
Further increasing s, we find the interference fringes be- 
come blurred and disappear eventually (see (a4)-(a6)), 
similar to the experimental observation in Ref. [25| . The 
gradual disappearance of the interference fringes can be 
partially explained by the increasing of the random rel- 
ative phase [25| during the SF-QC transition. In Figs. 
2(a5)-(a6), except for the density fluctuations along z 
direction associated with the random relative phase be- 
tween different subcondensates, we also see significant 
density fluctuations along y direction. The simultaneous 
existence of density fluctuations in these two directions is 
further demonstrated in Fig. 3 (c) and (d) for a typical 
case of s = 30.7. These experimental results suggest that 
in the QC region above the BKT transition, both the ran- 
dom relative phases and spatial phase fluctuations play 
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important roles in the expanded density distribution. 

To qualitatively analyze the density fluctuations, we 
consider the following density probability distribution 



P(n) = 



S{n- dn/2) -S{n + dn/2) 

Sn ■ S'total 



(1) 



where Sn is the width of a density interval, S'totai is the 
total area of the region occupied by the atomic gas, while 
S (n) is the area of the region where the density larger 
than n. The averaged density is Ug = 1/5'totai- For con- 
venience, we define a dimensionless density as uq = n/us. 
The dimensionless density probability distribution is then 
Poino) = nsP{n). 

The right column of Fig. 2 displays Pq ('^■o) calculated 
from the corresponding density distribution. In each im- 
age, the pixel region chosen for the calculation of Pq ('^•0) 
is that occupied by the cold atoms, as enclosed by the 
red dashed ellipse. In all our calculations of Pq (no)j 
driQ is just the horizontal spacing of discrete data points. 
Due to the optical noise, Pq (no) can be nonzero even for 
negative riQ. This can be understood from the atomic 
density formula based on the absorption imaging sig- 
nal. Assuming nofc is the minimum negative density for 
nonzero Pq (wq), optical noise concentrates in the region 
of no6 < ^0 < |nob| in each Pq (no) plot. Thus, ng > |nob| 
gives the effective region where the optical noise is neg- 
ligible and Pq (no) reflects the true density probability 
distribution. 

The blue solid lines in Fig. 2(bl)-(b6) represent the 
following exponential density probability distribution 
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It is obvious that, with increased s (and hence IJ j J\ 
Pq (no) has a tendency to Pq . In Fig. 2(b6), we see that 
Po (no) agrees well with Pg in the whole effective region. 
Despite of the extremely complex many-body state in the 
QC region above the BKT transition, our results show a 
simple universal behaviour. 

We have numerically simulated our experiments to ex- 
plain the exponential density probability distribution. 
When a released gas has experienced a free expansion 
over a time of t, the 2D density distribution is given by 
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Here (/3fej_ {x,y,t) (fkz {z,t) is the wave packet of the tribution in transverse (x-y) directions, and \ipkz {z,t) p 
expanded subcondensate initially in the fcth lattice site, the Gaussian density distribution along z direction, 
with \ipki_ {x,y,t) p being the Thomas- Fermi density dis- Two classes of phase fluctuations enter the expression. 



4>k± {x, y, t) represents the spatial phase fluctuations in 
transverse directions, while (pl {z,t) comprises the rela- 
tive phase fluctuations between different subcondensates. 
iVfc denotes the atomic number of the fcth subcondensate. 
The term 71,2^ is added to account for the optical noise. 
In numerical calculations, it was simulated with the re- 
alistic optical noise distribution in the region having no 
atoms. 
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FIG. 4: Theoretical simulation. Left column; simulated den- 
sity distributions based on Eq. ([3}. (al), (a2) and (a3) are 
three typical cases in which an atomic gas initially stays in 
SF, Type-I QC and Type-II QC regions, respectively. In (al), 
s = 11.2, S(l)k± = 0.27r. In (a2), s = 30.7, both 0^ and 0fc_L 
are completely random. In (a3), s = 30.7, 4)k± = 0, but 
(j)l is completely random. Each of the color map has a size of 
0.9x0.9 mm'^ . Right column: Black lines are the density prob- 
ability distribution calculated from the corresponding density 
distributions in the left column. Red lines are obtained un- 
der the same conditions except that the optical noise is not 
included. The exponential density probability distribution 
(blue lines) is also plotted for comparison. The data in (b2) 
clearly shows an exponential density probability distribution 
in the QC region above the BKT transition. 

Figs. 4(al) and (bl) give the simulated density distri- 
bution and density probability distribution for an atomic 
gas initially in the combined trap with s = 11.2. Spa- 
tial phase fluctuations {<pk±) have been assumed to be 
zero. Figs. 4(a2) and (b2) display similar calculations for 
s = 30.7, but with completely random phases 0fcj_ and 
0^.. The theoretical results in Figs. 4(b2) match the ex- 
ponential density probability distribution, in agreement 
with our experimental results. The coexistence of two 



random phases of (j)k± and 01 makes the system similar 
to the situation of the interference of randomly scattered 
waves |26[, where exponential distribution is also found. 
More general analytical derivations for Eq. ([3]) and the 
exponential density probability distribution are given in 
the Appendix. Figs. 4(a3) and (bS) give the simulation 
of s = 30.7 for a cold atomic gas in the QC region be- 
low the BKT transition. Accordingly, (/)fe_L is assumed to 
be zero, while there exists a completely random relative 
phase in (/)|. No exponential density probability distri- 
bution is found in this simulation. In addition, we do 
not notice regular interference fringes either. Note that 
in Ref. [23|, the QC region far below the BKT transition 
was studied experimentally. However, because there are 
only about 30 subcondensates, high-contrast interference 
fringes were still observed. Since our calculation is based 
on the experimental parameters in this work, the total 
number of subcondensates is much higher (up to 250). 

To further explain the experimental phenomena, we 
stress that there are two completely different situations 
in determining the BKT transition temperature for cold 
atoms in ID optical lattices. 

(1) BKT transition in QC and MI regions: In 
these two regions, the quasi-2D Bose gases in different 
lattice sites can be regarded as independent subconden- 
sates as far as the BKT transition is concerned. By min- 
imizing the free energy F = E — TS {E and S are the 
energy and entropy of the Bose gas in a lattice site when 
vortices are considered) , one can get the BKT transition 
temperature Tbkt = T2dI^ [2J|. In the MI region above 
the BKT transition, Eq. ([3]) can be used to simulate the 
density distribution by including a completely random 
relative phase in t/)^ [23| , if the averaged particle number 
per lattice site is much larger than 1. Thus, in the MI 
region above the BKT transition, one still expects the 
exponential density probability distribution for a large 
number of more than one hundred independent subcon- 
densates. 

(2) BKT transition in SF region: In the SF region, 
although the Bose gases in the lattice sites become quasi- 
2D, they are highly correlated. For a series of M highly 
correlated quasi-2D Bose gases, the free energy becomes 
F = ME — TS. If there are vortices in the quasi-2D sub- 
condensates, vortices in different lattice sites are highly 
correlated in spatial locations. Thus, S can be approxi- 
mated as the entropy of a single subcondensate. Then the 
BKT transition temperature becomes Tgj^j, — MTbkt- 
Since M > 100, we see that Tgj^j, >> T2d- This means 
that the temperatures of quasi-2D condensates are al- 
ways much lower than the BKT transition temperature. 
For the situation of s = 11.2 in our experiment, based 
on the estimation of the particle number fluctuations 
5Nk [131 and the relation of 5(1)1 - 2n/6Nk [2_5|], we de- 
duced a phase fluctuation of Scfif. ~ 0.2n « 2tt, which 
shows strong correlation between different lattice sites. 
As shown in Fig. 3(b), there are no noticeable density 



fluctuations along y direction, which confirms that the 
initially confined gas is well below the BKT transition 
even though T is already beyond Tbkt- 

We should also emphasize that, in the QC region above 
the BKT transition, the system cannot be regarded as 
a classical thermal cloud although the temperature is 
higher than the critical temperature T^. A classical ther- 
mal cloud features density fiuctuations as (5^n « n. In 
contrast, the exponential density probability distribu- 
tion corresponds to much larger density fiuctuations of 
(5^n « n^ [28|. This exponential density probability dis- 
tribution physically originates from the superposition of 
quantum sates with a completely random phase, while 
the universal behaviour lies in that the distribution is al- 
ways exponential for completely random superposition of 
quantum state [28| . The present work supports the long- 
standing belief that the quantum states in the QC region 
are strongly correlated, and there should be a simple uni- 
versal behaviour, once appropriate physical quantity is 
found 0^. 

H. W. thank the cooperation with Biao Wu about the 
dynamical universal behaviour of quantum chaos, which 
stimulates the present work. This work was supported by 
National Key Basic Research and Development Program 
of China under Grant No. 2006CB921406 and NSFC 
under Grant No. 10634060. 



[I] J. Hertz, Phys. Rev. B 14, 1165 (1976). 

[2] S. Sachdev, Quantum Phase Transitions (Cambridge 

Univ. Press, Cambridge, 1999). 
[3] P. Coleman, A. J. Schofield, Nature 433, 226 (2005). 
[4] P. Gegenwart, Q. Si, F. Steglich, Nature Phys. 4, 186 

(2008). 
[5] M. Imada, A. Fujimori, Y. Tokura, Rev. Mod. Phys. 70, 

1039 (1998). 
[6] S. Sachdev, Rev. Mod. Phys. 75, 913 (2003). 
[7] I. Bloch, J. Dalibard, W. Zwerger, Rev. Mod. Phys. 80, 

885 (2008). 
[8] M. P. A. Fisher, P. B. Weichman, G. Grinstein, D. S. 

Fisher, Phys. Rev. B 40, 546 (1989). 
[9] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, P. 

ZoUer, Phys. Rev. Lett. 81, 3108 (1998). 
[10] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, L 

Bloch, Nature 415, 39 (2002). 

[II] V. L. Berezinskii, Sov. Phys. JETP 34, 610 (1972). 

[12] J. M. Kosterlitz, D. J. Thouless, J. Phys. C 6, 1181 

(1973). 
[13] A. Polkovnikov, E. Altman, E. Demler, Proc. Natl Acad. 

Sci. USA 103, 6125 (2006). 
[14] Z. Hadzibabic, P. Kriiger, M. Cheneau, B. Battelier, J. 

Dalibard, Nature 441, 1118 (2006). 
[15] R. B. Diener, Q. Zhou, H. Zhai, T. L. Ho, Phys. Rev. 

Lett. 98, 180404 (2007); Q. Zhou and T. L. Ho, Phys. 

Rev. Lett. 105, 245702 (2010). 
[16] B. Capogrosso-Sansone, N. V. Prokof 'ev, B. V. Svis- 

tunov, Phys. Rev. B 75, 134302 (2007). 
[17] Y. Kato, Q. Zhou, N. Kawashima, N. Trivedi, Nature 



Phys. 4, 617 (2008). 
[18] N. Gemelke, X. B. Zhang, C. L. Hung, C. Chin, Nature 

460, 995 (2009). 
[19] K. R. A. Hazzard, E. J. Mueller, arXtv: 1006.0969 

(2010). 
[20] W. Kohn, Phys. Rev. 115, 809 (1959). 
[21] P. Pedri et al, Phys. Rev. Lett. 87, 220401 (2001). 
[22] S. Burger et al, Europhys. Lett. 57, 1 (2002). 
[23] Z. Hadzibabic et al, Phys. Rev. Lett. 93, 180403 (2004). 
[24] C. J. Pethick, H. Smith, Bose-Einstein condensation in 

dilute gases. (Cambridge University Press, Cambridge, 

2008). 
[25] C. Orzel, A. K. Tuchman, M. L. Fenselau, M. Yasuda, 

M. A. Kasevich, Science 291, 2386 (2001). 
[26] P. Sheng, Introduction to Wave Scattering, Localization, 

and Mesoscopic Phenomena. Chap. 10 (Academic Press. 

Inc., California, 1995). 
[27] J. Javanainen, M. Y. Ivanov, Phys. Rev. A 60, 2351 

(1999). 
[28] H. Xiong, B. Wu, arXiv: 1007.2771 (2010). 

Appendix 

Bose-Hubbard model. The combined potential of a 
one-dimensional (ID) optical lattice and a harmonic trap 
is given by 



V — —mu I (x 



y^) + -ttujjIz^ + sEr sin^ i -^ 



(4a) 
Here Er is the recoil energy of an atom absorbing one lat- 
tice photon with wavelength A. For an atom in a lattice 
site, it experiences an effective harmonic potential along 
z direction with angular frequency uiz — 2^Eji/h. 

Bose-condensed gases in ID, 2D and 3D optical lattices 
are widely studied by the following Bose-Hubbard model 
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(4b) 

Here rii — afii is the number operator at the ith site, 
with aj (fli) being the boson creation (annihilation) op- 
erator. U and J represent the strengths of the on-site 
repulsion and of the nearest neighbor hopping, respec- 
tively. Si describes an energy offset due to the harmonic 
trap, and /j, is the chemical potential. The phase diagram 
given by Fig. 1 in the text is calculated with the above 
combined potential and Bose-Hubbard model. 
Universal exponential density probability distri- 
bution. The many-body quantum state in the QC region 
can be written as in a general way 



(4c) 
Here Nk denotes the atomic number in the fcth lattice 



site. The total atomic number is then N 



^k—k^ 



Nk 



with 2kM + 1 being the total number of occupied lat- 
tice sites. For strongly correlated quantum gases in 
the QC region, the function C{Ni,- ■ •,A^2fcM+i) is ex- 
tremely complex, and no analytical expression exists. Af- 
ter switching off the combined potential, the density dis- 
tribution is then 

n {x, y, z, t) - (*,,, {t)\ $t$ 1*^^^ (t)) . (4d) 

Assuming that the initial wave function 
in the fcth lattice site is ipk {x^y,z,t = 0) = 
'/5fe_L (x, y,t = 0) ipkz [z, t = 0), for long-time evolu- 
tion so that the expanded subcondensates have sufficient 
overlapping with each other, we have 



n{x,y,z,t) 



^v^e**^(^'*Vfc±(a;,y,0 



\VOziz,t)\' 



(4e) 

This formula holds in the QC region, because there 
are non-negligible particle number fluctuations in each 
lattice site. Here 4'1{z,t) = 0^ + 0°, with 0^ being 
a completely random relative phase in the QC region. 
(/>" = ra{z — kX/2) /2ht can be obtained directly from 
the free expansion along z direction. Without the ran- 
dom phase (j)^. , the phase ip^ will lead to clear interference 
fringes along z direction with a period of Anht/m\ 2\ . It 
is understood that there would be no interference fringes 
along z direction if 0| is a completely random relative 
phase, as shown in Figs. 2(a5)-(a6) in the text. 

In the QC region above the BKT transition, there are 
two sorts of random phases: (i) the completely random 
phase <j)\ (z, t) in the QC region; (ii) the spatially relevant 
random phase (j)ki_ {x, y, t) above the BKT transition. In 
this region, the density distribution can be further writ- 
ten as 



n{x,y,z,t) 



Ev^ 



^ict>k±ix,y,t) i(t>l{z,t) 



\vo±ix,y,t)(poz {z,t)\ 



(4f) 



Here (/)fc± {x,y,t) physically originates from the spatial 
phase fluctuations of quasi-2D Bose gas above the BKT 
transition. 

The density distribution recorded by a CCD is a 2D 
density distribution n2D iy,z,t) — J dxn{x,y, z,t). By 
assuming that \ipQ± (x, y, t) (poz (z, t)\ is locally constant 
(i.e. the local density approximation), after integrating 
out the X variable, we have 
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n2Diy,z,t) 

a/ V.tCi, t -+- I ' 

k 



V2k 



L=Y,c^k{y,z,t)e^^''(y'^'''> 
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iy,z,t), 



(4g 



with n%{y,z,t) = J dx VNipo± {x,y,t) (poz {z,t) . 
ak{y,z,t) and /3k{y,z,t) are statistically independent 



random real function, because they originates from 
two different random phases 4)l.{z,t) and (f>k± {x,y^t). 
Pk {y, z, t) distributes uniformly between — tt and tt. The 
normalization condition requests (a^) = 1- 

Based on the central limit theorem [6[, we get the fol- 
lowing exponential density probability distribution 

g-n2D(x,y,t)/n°^{x,y,t) 

P{n2D{x,y,t)) = ^- . (4h) 

nlD{x,y,t) 

For uniform n^jj, it is straightforward to get Eq. (2) 



in the text. For nonuniform 



'2_D' 



with the local den- 



sity approximation, Eq. (2) in the text is a very good 
approximation obtained from Eq. (j4h|) . 

In Figs. 2(b5)-(b6) in the text, we see that the density 
probability distribution agrees well with the exponen- 
tial density probability distribution. However, because 
the system temperature is slightly smaller than the crit- 
ical temperature for Fig. 2(b5) (s = 22.3), we expect 
Fig. 2(b6) agrees better with the exponential density 
probability distribution. Here we give the semilog plot 
of the density probability distribution for s = 22.3 and 
s = 30.7, respectively. We do find that the situation 
of s = 30.7 agrees better with the exponential density 
probability distribution. 
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FIG. 5: Semilog plot of the density probability distribution 
for s = 22.3 and s = 30.7. The blue line gives the exponential 
density probability distribution. 

System temperature. In our experiment, the total 
atomic number can be accurately measured. The atomic 
numbers of Figs. 2(a4)-(a6) in the text are 1.10 x 10^, 
1.19 X 10^ and 1.13 x 10^, respectively. To drive the 
system into the deeper region of superfluidity, we have 
to further lower the temperature of an atomic sample 
before loading it to the lattice. This is implemented by 
kicking more atoms out of the magnetic trap during the 
evaporative cooling stage. Specifically, for Figs. 2(a2) 
and (a3) in the text, the atomic numbers are reduced by 
about a factor of two compared with Figs. 2(a4)-(a6) in 
the text. 

However, the measurement of the system temperature 
is hampered by the lack of accurate thermometry in opti- 



cal lattices. Actually, the ratio T/Tc is a more important 
quantity than T in characterizing the transition from su- 
perfluid to quantum critical regions. Therefore, the inac- 
curacy of system temperature does not affect the physics 
of the present work, as long as T/Tc can be somehow 
determined with a satisfying accuracy. 

As predicted in Ref. [3| and demonstrated in a recent 
experiment [^ , the appearance of superfluid in an inho- 
mogeneous lattice is associated with a bimodal structure 
in the density profile. The condensate fraction N^/N 
can be obtained by fits to the bimodal distribution. As 
we have known, Nq/N shows a temperature dependence 
having a characteristic shape as 



No/N=l-iT/T,r. 



(4i) 



Apparently, T/Tc can be directly determined if the pa- 
rameter a is known. Blakie et al. [5| had performed 
numerical calculations of the critical temperature as well 
as the condensate fraction, for a combined trap similar to 
that in our work. Their results show that a is close to 3 at 
a moderate lattice depth, but it may drop slightly with 
increasing lattice depth. Nevertheless, at a high depth 
level, the condensate fraction is usually very small, and, 
as a result, the value of T/T^ derived from Eq. (|4i| be- 
comes insensitive to the potential errors of a. Therefore, 
we just set a = 3 in the calculation of T /Tc, which seems 
to be a reasonable assumption. 



FIG. 6: A bimodal fit to the interference pattern in Fig. 
2(a3) in the text. The total atom number is measured as 
A^ ~ 5.9 X 10*, and the atomic cloud was released from the 
combined trap with a lattice strength of s = 11.2. Black line, 
the measured atomic density which is obtained by integrating 
over each column of pixels. Red hne, a fitting curve including 
six Gaussian peaks, where the three narrow peaks at the up- 
per part are attributed to superfluid. The condensate fraction 
is A''o/iV « 0.58, corresponding to T/Tc ~ 0.75. 



For each absorption image, we first obtain the atomic 
density along the z dimension by integrating pixels in 
each column. The narrow peaks riding over broad ones 
are regarded as the condensate parts. We then fit both 
the broad and the narrow peaks using Gaussian profiles, 
so as to get the condensate fraction. Apparently, this 
method works only when T < T^- In Fig. 2(a6) in the 
text, the interference peaks vanish completely, indicating 
that T > Tc- In this case, we are unable to extract the 
quantity T/Tc any more with this method. The universal 
exponential density probability distribution found in this 
work also implies that the accurate measurement of sys- 
tem temperature in the QC region above BKT transition 
is difficult. 
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